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Abstract 

We use our model of the small scale structure on cosmic strings to develop further 
the result of Siemens, Olum, and Vilenkin that the gravitational radiation length scale 
on cosmic strings is smaller than the previously assumed TGut. We discuss some of the 
properties of cosmic string loops at this cutoff scale, and we argue that recent network 
simulations point to two populations of cosmic string loops, one near the horizon scale 
and one near the gravitational radiation cutoff. 



1 Introduction 



Cosmic strings and their dynamics have been studied for more than twenty-five years. Al- 
though much is known about the formation and evolution of cosmic string networks [TJ [2] , 
the complexity of these systems still stands in the way of a complete understanding. It 
is generally believed that, after their formation, the networks evolve into a scaling regime 
where all length scales grow linearly with the Hubble time. This would imply that the typical 
size of cosmic string loops is a fixed fraction a of the horizon scale]]] However, there is no 
consensus on the value of the proportionality constant a; in fact, estimates range over tens 
of orders of magnitude. 

A key question is whether the smoothing of the string by gravitational radiation is nec- 
essary for scaling of the loop sizes. If so, a will depend on the dimensionless string tension 
Gfi, because the emission of gravitational radiation is proportional to this parameter. If not, 
then the evolution of the network will be purely geometric, and a will be a pure number, 
independent of Gfi. For many years, simulations appeared to show loops forming at the short 
distance cutoff scale, and this was interpreted as implying the need to include gravitational 
radiation. Several recent simulations, however, indicate that the loops are actually forming 
above the cutoff scale, so that a will be independent of Gfi. Refs. [H E] finds notably large 
loops, a ~ 0.1, while Refs. [6j [JJ find loops a few orders of magnitude smaller. 

In Ref. [8] we attempted to construct an analytic model of the small scales in cosmic 
string networks. We had some success in understanding the two-point functions of the 
tangent vector and velocity on long strings. However, loop production was found to diverge 
at small scales, and we had no means to deal with the resulting strongly nonlinear dynamics. 

In Sec. 2 we use our model to study gravitational radiation from long strings. It was long 
assumed that this radiation smoothed the strings on scales of order the horizon length times 
Gfi, so this would be the size of loops if they formed at the gravitational radiation cutoff [9j 
ITU] . However, Ref. [T3] showed the radiation from long strings had been overestimated, 
and so it becomes important only at a shorter length scale, proportional to a larger power of 
Gfi [E] ; the exponent depends on the power spectrum of fluctuations on the long string. Here 
we rederive the results of Ref. [TBI [H] in perhaps a simpler way. Our result for the exponent 
of Gfi differs somewhat from Ref. [2] . This difference arises in going from the periodic wave 
train calculation of Ref. [J3] to the random distribution on the actual network. Also, we use 
input from our model [8] and from simulations [7] to determine the actual power spectrum 
on the long string. 

In Sec. 3 we discuss loop production. We use the results from Sec. 2 to find the cutoff 
on loop sizes. It is notable that in our model these smallest loops will be highly relativistic. 

1 There is not complete agreement even on this point: see Ref. [3] and references therein. 
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We then consider the recent simulations of Ref . jU [5] , which are able to cover the longest 
time scales, and which show large loops. These simulations also show a population of small 
loops near the cutoff scale, which are claimed to be transient. We instead argue on several 
grounds that these smaller loops are a real effect, so that the true loop distribution has two 
peaks. One scales at a ~ 0.1, while the second is presumably at the scale set by gravitational 
radiation. We discuss how the resulting picture affects current limits and future sensitivity 
of searches for cosmic strings. In particular, Refs. [I5J [TBI [T7] showed that the a ~ 0.1 loops 
produce a much larger stochastic background than in previous network scenarios, tightening 
existing bounds and increasing the potential for discovery. If these loops represent only a 
fraction of the total string length then this effect is somewhat reduced, but still significant. 
We also discuss the gravitational wave signatures of the small loops, and the interaction 
between the short distance structure and the string cusps. 

2 Gravitational radiation 

2.1 Model for small scale structure 

Let us first review some of the results of Ref. [H] for the small scale structure on a long cosmic 
string. The motion of a cosmic string in a three-dimensional target space can be described 
by a map x(cr, t), where a parametrizes the distance along the string and t is the FRW time. 
Denote d a by a prime and a(t)dt = d T by a dot (where a(t) is the scale factor and r is the 
conformal time), and define e = a/x /2 / (1 — x 2 ). Then in transverse gauge, x • x' = 0, the 
quantities p± = x ± -x' are unit vectors. In flat space these describe simply left- and right- 
moving waves. In an expanding universe the equations of motion derived from the Nambu 
action introduce a correction to the wave equation which is proportional to a/ a. 

To better understand the small scale behavior of cosmic string networks a model was 
developed in [8] that focused on the two-point functions (p±(cr, t) ■ p±(o"', t)). This correlator 
must vanish at scales larger than the horizon size due to causality, while it encodes the short 
distance structure observed at smaller scales. Assuming that stretching of the network as 
spacetime expands gives the dominant behavior at scales small compared to the horizon but 
long compared to the gravitational radiation scale, the two-point function was determined 
to be 

(p ± M)-p±(^*)>~i-^A) 2x , (2-1) 

where I = a(t)e(t) \a — a'\ is the physical separation (weighted by the local value of 7 along 
the string). This quantity is a function of the ratio l/t only and so it scales properly. The 
exponent x is a function of the rate of change of the scale factor and also the RMS velocity 
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of the strings, v 2 . For a power law expansion of the form a(t) oc t u one finds [8] 



va. 



X = ; , a = l-2v 2 . (2.2) 



In a flat spacetime the virial theorem imposes v 2 = 1/2 (in natural units) but due to 
expansion v 2 < 1/2. Thus, for a decelerating universe the exponent x is a positive number, 
taking the specific values x ~ 0.10 in the radiation dominated era and x ~ 0.25 in the matter 
dominated era. 

These exponents are in reasonable agreement with the two-point functions given in Ref. [7] 
over a range of scales below the horizon size, with the normalization A found to be ~ 0.60 
in both the matter and radiation epochs. At smaller scales the simulations seem to indicate 
a larger exponent x ~ 0.5, corresponding to the functions p± mapping out a random walk 
on the unit sphere. We are unable to identify a physical mechanism that would produce 
this effect, and believe that the discrepancy is due to transients in the simulations. In any 



event, our discussion below just uses the general power law form (2.1), and one can insert 
any assumed values for x an d A. 

Also useful for the computations in the next section is the small fluctuation approximation 
introduced in [8] . In this framework the structure on a small segment is expressed as a large 
term that is constant along the segment plus a small fluctuation: 

p + (a,t)=P + (t)+w + (a,t)-^P + (t)w 2 + (a,t) + ... , (2.3) 

with P+{t) 2 = 1 and P+(t) • w + (cr, t) = 0. An equivalent equation can also be written for the 



right-moving vector p_. Equation (2.1) can then be expressed in terms of the fluctuations 

as 

1 - (P±M) ■ P±(</,t)> = l([w ± (a,t) - w ± K,t)] 2 ) + 0(wi) « A{l/t) 2 * . (2.4) 
2.2 Radiation and back-reaction 

Gravitational radiation is important on scales short compared to the Hubble time and so we 
can use the flat metric diag(— 1, +1, +1, +1), and also make the gauge choice at = lj^] The 
modes are then functions only of u = t + o or v = t — a, so that 

(p+(«) • p+K)> ~ 1 - A[{u - u>)/t} 2 * , <p_( V ) • p_(«0> « 1 - A[(v - v')/t] 2 * ; (2.5) 

(the explicit t on the right-hand side is effectively a fixed parameter, varying only over the 
Hubble time). 



2 This choice is possible when we consider scales small compared to the Hubble time. In Sec. 2.4, when 
we again consider cosmological evolution, we must reintroduce a(t) and e(t). 
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The energy radiated per unit solid angle in the k direction is given by [TU H21 [13] 
dAE Gfi 2 r °° 



^ [™ duu 2 {\A\ 2 \B\ 2 + \A*-B\ 2 -\A-B\ 2 } , (2.6) 
3tt J 



dQ 16? 

with the left- and right-moving contributions given respectively by 

/+oo ( i f u 

dup^{u) exp< — / du'k- p+(u') 
oo I 2 Jo 

"+0O 



dv f_ (v) exp< — / dv'k ■ V- 
oo I 2 ,/ 



(2.7) 



Here, is the gravitational radiation wave vector, and p± are null 4-vectors whose time 
component is identically 1 and with p± being the unit vectors above. These are related to 
the standard notation 

x»(u,v) = l -[a»{u) + W{v)} (2.8) 

via p%(u) = a^'(u) and p1(v) = b^(y). In the spirit of the small fluctuation approximation 
we orient a given segment of string mainly along the z-axis and consider fluctuations w± in 
the perpendicular x?/-plane. Thus, 



ft(u)~ ll,W + (u),yJl-W 2 + (u)\ , P 1(V)~ ll,W_(v),-yJl-wl(v)\ . (2.9) 



We can ignore the last two terms in equation (2.6): they are equal when B^ is real, and so 



cancel when ensemble-averaged over the short distance structure. 

Let us review the argument of Refs. [TBI HI]- First linearize the modes (2.7) in the 
oscillations, so the exponential factors become e lkaU and e lkbV respectively, with 

u = k a + k b , k z = k b -k a . (2.10) 

Here u and k z are the frequency and wavenumber of the radiation. Then the ensemble 
averages in the linearized approximation («) are 

/+oo 
due ik a u-e\u\ ((p+(M) . p+(0)) _ 1)= kAcx fc -(l+5*) f -5* , 
■oo 

(\B(k)\ 2 ) w i v Ac x k; (1+2x) r 2 x . (2.ii) 

Here c x = 2 sin(7rx) T(l + 2%), taking the values 1.25 and 0.57 in the matter- and radiation- 
dominated eras respectively. The convergence factor e~ £ '"' accounts for the decay of the 
correlations on horizon scales; its detailed form is unimportant at shorter scales. The factors 
l u and l v are volume regulators: we artificially cut the oscillations off to produce finite trains, 
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but these regulators of course drop out when we consider rates per unit time and length. 
Using dud cos 6 = 2dk a dk b / \k a + k b ) , we can write the total energy radiated as 

f..2 roc poo 

(AE) = dk a dk b (k a + k b )(\A(k a )\ 2 )(\B(k b )\ 2 ) 



4tt 

Gy 2 A 2 c 2 r°° [°° k a + k b 

"—tr'^Jo dka L dh • (2 ' 12) 

The volume of the world-sheet is V = ^l u l v , so this translates into a power per unit length 

dP\ Gn 2 A 2 c 2 x r°° i7 r°° k a + k 



i i dk„ \ dk b — — , „ — — . (2.13) 

dz I 2vr J a J (k a k b y+ 2 *t*x 1 ; 

The total energy of the wavetrains is 

^\ A (ka)\ 2 + »J Q ^\B(h)\ 2 - (2-14) 
Isolating a momentum range dk a , we have 

A\A(k a )\ 2 ^ -^\A(k a )\ 2 I™ dk b (\B{k b )\ 2 ) . (2.15) 
^ Jo 

We have used the fact that momentum conservation determines that the energy coming from 
the A and B modes is in the proportion k a to k b . A given point u interacts with the v train 
for a time |At), so the rate of decay becomes 

j t \A(k a )\ 2 « -G^k a \A(k a )\ 2 Ac x ^ ^i^- (2.16) 

The integral is dominated by the lower limit: most of the energy loss from the left-moving 
mode k a comes from its interaction with right-moving modes of much longer wavelength. 
Cutting off the lower end of the k b integral at the horizon scale, k b ~ we find that the 
decay rate is of order Gfik a . This is faster than the Hubble time for 

k a > 0(G/it)" 1 , (2.17) 

meaning that modes with wavelengths A < 0(G/j,t) are exponentially suppressed. Thus, 
Eq. (2.17) reproduces the usual (naive) estimate for the gravitational damping length. 



2.3 Radiation and back-reaction improved 

The calculation above implies that the energy loss from the left-moving modes k a comes 
primarily from their interaction with much longer right-moving modes at the horizon scale. 
Ref. |T3] argued that this was paradoxical, because a wave that encounters an oncoming 
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perturbation with much longer wavelength is essentially traveling on a straight string, and 
such a wave will not radiate. They showed that this paradox arose due to the neglect of the 
exponential terms in the modes (2.7). With their inclusion, the radiation is exponentially 
suppressed when the ratio of k a to k b becomes sufficiently large. 
Ref. 1 13 1 considered monochromatic wavetrains 



w+(u) = e a cos(k a u) , w x _(v) = e b cos(k b v) , (2.18) 
and found that the gravitational radiation was suppressed unless 

k b /e 2 b >k a , h<k a /e 2 a . (2.19) 



The first of these relations will cut off the integral (2.16) below the horizon scale, but first 
we need to extend it to the incoherent spectrum on a long cosmic string. We will do this in 
a systematic way below, but we can anticipate the answer. For a continuous spectrum, of 
course, a single frequency makes a contribution of measure zero. The natural extrapolation 
is then to replace e a , e b with the RMS average in bands of frequencies k ~ k a , k b . That is, 

4- r d ^^(k)-*-(kr 



k b /A 2VT Jk b /A 2VT 



CO 



/\i2 



1 

~2 

(hi)-' 2 * ; (2.20) 



dv I dv' g(v)g(v')*[w-(v) — w_(t/)] 

oo 

2X . 



here A = ^fe is a plausible choice, corresponding to integrating over one unit of log frequency, 
but the precise value does not matter. In the second line the filter function 

r Akb dk 

g(v) = / %-e*° (2.21) 
Jk b /A 2vr 

has a width of order l/k b and a magnitude of order k b . In the third line, we have used the 
two-point function (2.4). Thus, 

{k b t) 1+2x > k a t . (2.22) 



Using this lower cutoff in the decay rate (2.16) gives 

±\A{k a )\ 2 ~ -Gn{k a t)V^h-M{K)\ 2 - • (2.23) 

dt r GR (k a ) 

This is faster than the Hubble rate for 

K^G^-^H- 1 . (2.24) 

The gravitational length scale is reduced from the naive 0(Gfit) by a factor (Gfi) 2x . The 
qualitative conclusion is the same as in Ref. [14J, but the exponents do not seem to agree. 
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The result there was k a > (Gfi) ( 1+2/3 )/ 2 t 1 . It is not clear whether (3 is to be identified with 



X + \ (as suggested by Eq. (21) of Ref. [14J, compared with Eq. (2.16) above) or with x 



(as suggested by [H] Eq. (31)), but in either case the exponents differ. This stems from a 
different method of converting from the single-mode result to the continuous spectrum; our 
exponent will be borne out by the more formal treatment below. 

Finally, recalling Sec. 2.1, our analytic model j8] gives for the exponent 1 + 2% the values 
1.2 in the radiation era and 1.5 in the matter era; using instead the simulations [7] would 
yield an exponent around 2.0 at the shortest scales in both erasj^] 

The necessary correction to the naive radiation formula comes from the previously ne- 



glected exponential factor in B fM (k), Eq. (2.7) [13]. Expanding the exponential to second 
order in the fluctuations gives 



\B{k)\< 



/+oo f+OO 
dv / dv' (p-(u) • P-(f') — 1) 
-oo J — oo 



exp < ik h {v' — v ) — 



ik i 



dv"w 2 _(v") } . (2.25) 



Noting that k\ = 4k a kb, when the first term in the exponent is of order one, the second and 
third terms in the exponent are respectively of order (k a t) 1 ^ 2 (k b t)~ 1 ^ 2 ~ x and k a t{k b t)~ l ~ 2x . 



Thus, as kb decreases they become important at precisely the scale (2.22) where the small 



amplitude calculation is expected to break down. The linearized form for A(k) remains valid, 



so we can simply insert the corrected form (2.25) into the decay rate (2.15) 



By completing the square, the exponent in Eq. (2.25) can be written as 

. .2 

i — I ki , a I dv"w^(v") 



V 



ik 2 



[V — V 



r dv" f V dv'" [w_(u") - ™-{v'")f . (2.26) 

J V J V 



Converting J dk^ — > (4/c a ) 1 J dk\ — > (4:irk a ) 1 J d 2 kj_, we have at fixed k a the integral 

7 i(p_(«). P _(«')-i) 



dk b \B(k)[ 



dv 



dv' 



V — v 



%k? 



exp 



K 1 "' - v ) Jv 



dv" 



dv"'[w_(v")-w4v"')} 2 } . (2.27) 



We have used — k z m k a ^> k^, since this is the regime where the correction is important. 



3 Numerically, there are two changes from Ref. |14j : the corrected expression for the exponent, and a more 
accurate estimate of the power spectrum. Ref. [2] effectively estimated the latter using v 2 — 0, so that their 
Eq. (31) top = y\* =a = v/{l-u). 
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We now approximate, replacing both the exponent and the prefactor with their mean 
values from Sec. 2.1. Then 



dk b (\B(k)\ 2 ) = 2M v lm 



+00 V 2 X -1 

dv 



2l v b x (k a t) * ^~ . 



exp 



ik a A 
4v 



dv" / dv 



\v" -v"'\ 2 x 



t 2 * 



Here 



y = 2 sin — 

l l + 2 X 



r 



2X 



[4(1 + x )] 2 x/d+2x) 



A 



V(l+2x) 



(2.28) 



(2.29) 



is 2.5 in the matter era and 2.1 in the radiation era. We can improve this approximation if 
we make the assumption that the ensemble is gaussian. The first order correction, keeping 
contractions between the prefactor and the exponent, involves a straightforward calculation; 
its inclusion simply renormalizes the constant b x — > (1 + € x )b x , with 



r(i + 2 X ) 2 1 

T(2 + 4 X ) 2(1 + 4 X ) 



(2.30) 



However, this represents a small correction since e x = 0.035 in the radiation epoch and 



0.045 in the matter epoch. Thus, to good approximation, 



1 



TGR(ka) 



(2.31) 



This is faster than the Hubble rate for the scales (2.24), as deduced earlier 



2.4 The long-string two-point function 

We can now improve our earlier result for the two-point function [8] through the inclusion 
of gravitational radiation. Define 



H(k, t) = a(t)e(t) / da e~ iKa (w+(a, t) ■ w + (0, t)) , 



(2.32) 



where a is related to the actual length along the string by I = a(t)e(t)a; here a = t u and 
e(t) = a~ 2 '" 2 . Combining the stretching given by Eq. 2.26 of Ref. [8] with the gravitational 
radiation found above gives 



^-H(n,t) = --aH(K.t) / . . H(k, t) 

dt V ; a V ; T GK {K/ae) V ; 



(2.33) 



This is readily integrated, to give 

\nH( K ,t) = -alna(t) + ln/(/«) - b x (l + x)(l + 2x)G/iK 1/(2x+1) t 1/(1+x)(1+2x) . (2.34) 
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log (l/t) 



Figure 1: The effect of the inclusion of gravitational radiation on the long-string two-point 
function is to smooth the string on scales below ~ 20(Gfi) 1+2x t. The solid line corresponds 
to the matter epoch and the dotted line to the radiation epoch. The value Gfi = 10~ 9 has 
been adopted for illustrative purposes. 



The function f(n) is obtained by matching onto the known result (2.4) at early times, 
giving f(n) = Ac x k~ 1 ~ 2x . Expressing the result in terms of the physical momentum scale 
k = n/a(t)e(t) gives 

J{k,t) = H(n,t) = Ac^ikty^exp [-b x {l +x)(l + 2 X )Gfi{kt) 1/{1+2x) ] . (2.35) 

Thus, unless we consider modes with such short wavelength that the exponential sup- 
pression factor kicks in, we have J(k) ~ k~^ l+2x \ At length scales smaller than ~ (Gfj>) 1+2x t 
this behavior is altered as a result of the inclusion of gravitational radiation, in such a way to 
make the two-point function smoother. In Fig. [T] we show the numerical Fourier transform 
of J(k,t), in terms of the correlation p3 [8] 

^ l-( P+ (/,t)- P+ (0,t)) _ (u; + (/,t) 2 )-(w + (/,t)-w + (0,t)) 
1 C ° rr *~ 2(1 -v 2 ) ' 2(1 -v 2 ) • {I - 6b) 



The smoothing at short distance seen in Fig. [T] resembles that found in the numerical 
simulations of Ref. [7j, but here it is due to an additional physical effect, gravitational radi- 
ation, whereas in Ref. j7] such smoothing was found even without gravitational radiation]^] 
Again, it remains to be seen whether the result of Ref. [7j is a real effect or a transient. In 

4 As an aside, the correlator in Fig.[l]goes as I 2 at short distance, corresponding to an ensemble of smooth 
functions, whereas it appears to be close to I 1 in Ref. [7], corresponding to a random walk in p±. 
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Fig. [T]we have assumed the latter, so that the two-point function follows a simple power law 
down to the gravitational radiation scale. 

3 Loops 

3.1 Small loops 

The exponential falloff in the world-sheet two-point function at large wavenumber implies 
that the strings are very smooth at short distance. Inspection of Fig. [T] indicates that this 
smoothing sets in at 

l c « 20(G^) 1+2x t (3.1) 

in both the matter and radiation epochs. The exponent we obtain is smaller than the 
estimates in [Hj: 1.2 in the radiation dominated era and 1.5 in the matter dominated era 
here, compared to 3/2 and 5/2 respectively. The difference is primarily from using a more 
accurate model of the effect of stretching, leading to a power spectrum with a slower falloff. 
The length l c is the approximate size of the smallest loops. In particular, the divergent loop 
production found in Ref. [S] is cut off at this scale. 

It is customary to parameterize the size of small loops as I = eTGfit with T ~ 50. The 
parameter e corresponds to the lifetime of the loop in units of the FRW time, if gravitational 
radiation is the principal decay mode. Thus we have found 

e « 0.4(G» - 5 (3.2) 

in the matter era, while for the radiation era 

e w 0.4(G» ' 2 . (3.3) 

However, for the phenomenology of the small loops there is another important property, 
first noted in Ref. [8]: they move extremely fast. The point is that the functions p+(u) and 
P-(i>) must have equal mean values on the loop (the condition for the loop to form). If the 
loop is small, so these functions have limited range, this implies that both remain near one 
point on the sphere, and the velocity (the average of p+(u) and P-(i>) over the loop) will be 
very close to unity. From Eq.(4.34) of Ref. [8], the Lorentz contraction factor is 

7 = (//t)^[(l + X )(l + 2x)M] 1/2 . (3.4) 

This is ~ 0.8(Gyu)~ 0,37 in the matter era, and is of order 10 3 for tensions of interest; in the 
radiation era it is ~ 1.0(G/i)~ 012 and is of order 10 1 . Since this factor is so large, we must 
be careful to recall that the specific definition of length is 

dl = a(t)e{t)da = a(t) ^'7(1 - ± 2 )da (3.5) 
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and so I is 7 times the rest frame length Z rest . It is equal to twice the inverse period of the 
loop in the FRW rest frame. 

To see one effect of the boost, in Ref. [17] the total stochastic gravitational wave spectrum 
from small loops was considered. The frequency spectrum simply scales inversely as the size 
of the loops if they are at rest, u> oc 1/Z res t- In the FRW rest coordinate this frequency is 
reduced by a factor of 7. However, the large boost pushes essentially all of the radiation into 
the forward direction, where the frequency is enhanced by a factor (1 — v )~ x ~ 2 , ~f 2 . Thus 
the gravitational radiation spectrum corresponds to an effective loop size 

e cff « 7" 1 e re st/2 « -y~ 2 e/2 . (3.6) 

This is e e fr ~ 0.3(G/i) L25 in the matter era, meaning that Z e ff ~ 15(Gyu) 2,25 . In the radia- 
tion era e e ff ~ 0.2(G/i) a44 , meaning that l eS « 10(G/z) L44 . We should note that Ref. [T7] 
considered the sensitivity in the whole (G/i, e) plane, and we are simply highlighting one 
line in this plane. It is notable that Advanced LIGO is more sensitive than LISA to these 
smallest loops, because their small size puts them above the LISA frequency range. In fact, 
for interesting values of the tension the loops will produce a periodic signal in the LIGO 
frequency band, which may be observable at Advanced LIGO [T8] . 

3.2 Loops: large and small 

Finally, let us discuss how these small loops might fit into the overall picture of the string 
network. In Ref. [8] we found that loop production from long strings diverges (as measured by 
total string length) at small sizes. To match the scaling value for the rate of long string loss to 
loops we would need to cut this divergence off at I ~ 0.18i; this is small in the sense of being 
an order of magnitude below the horizon length, but very large compared to the scales that 
we have discussed thus far. However, the emission of loops from long strings does not end the 
story, since the divergent loop production continues on the large loops themselves. We were 
unable to go further analytically, but conjectured that loops would fragment extensively until 
they reached a non-self-intersecting configuration at a scale we guessed to be a few orders of 
magnitude smaller, I ~ 10 3 t to 10 _4 t. This appeared to be consistent with the simulations 
of Refs. [61 [JJ . If this were the case then the gravitational radiation smoothing that we have 
discussed would be largely irrelevant, because loop production would already be cut off at 
much larger scales. 

However, Refs. [U E] have recently been able to simulate over longer time scales, and 
argue for a rather different picture. Namely, they find two populations of loops, one of which 
scales at a = l/t ~ 0.1, and a second which remains at the simulation cutoff length. They 
argue that a ~ 0.1 is the actual loop size, and that the small loops are a transient effect in 
the simulations and will not be present in the real networks. 
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We would argue that this interpretation of the small loops as transients is implausible. 
First, there is no sign in Refs. jH E] that the total fraction of string length going into the 
small loops is diminishing in time. Second, if this were a transient effect then its natural 
time scale should be small in proportion to the size of the loops. Instead, these loops are 
still being produced after several Hubble times, showing that some longer-distance process 
is continuing to feed the small loop production. Indeed, we have identified this mechanism 
already [S] : the small scale structure on long strings and large loops, which originates at the 
horizon length and then is carried down to smaller relative lengths by the expansion of the 
universe, leads to production of loops on arbitrarily small scales. 

Thus we are led to suggest the following picture. After a large loop is formed, production 
of smaller loops continues near its cusp regions. In the end, there is a large loop without 
cusps or self-intersections, reduced in size by the excision of the cusp regions. There is 
also a population of much smaller loops; based on the non-scaling seen thus far jH E], we 
conjecture that these will be at the gravitational radiation length]^ Inspection of Fig. 4 of 
Ref. [5] suggests that as much as 80% of the total string length goes into the small loops, 
with 20% remaining in the large loop. 

As has been discussed in Refs. [151 EE 02] > the a ~ 0.1 loops produce a much larger 
stochastic background than in previous network scenarios, tightening existing bounds and 
increasing the potential for discovery. If these loops represent only a fraction of the total 
string length then this effect is reduced. For loops produced and decaying during the radi- 
ation era, the total energy density today is proportional to (Gfxa) 1 ^ 2 . Thus, if only 20% of 
the string length goes into large loops, the bound on Gfi is weaker by a factor of 25. For 
example, the pulsar timing bound G[i < 8 x 10~ 9 for reconnection probability p — 1, shown 
in Fig. 2 of Ref. [17], would be weakened to Gji < 2 x 10~ 7 , very close to the current CMB 
bound [19]. Foreseeable improvements in pulsar timing, e.g. [20], will improve the sensitivity 
by several orders of magnitude in G[i. 

3.3 Cusps 

Our work also allows us to address an old question, the interplay between the small scale 
structure and the cusps [T5] . A cusp arises when the functions p+(u) and P-(f) cross on the 
unit sphere [2T] . The size of the cusp is the inverse of the 'velocity' of the functions p± when 
they cross, their rate of change with respect to the length along the string. Now, the average 
velocity of these functions on the long strings is (going again to the gauge e(t) = a(t) = 1) 

V 2 = (<9 CT p + (a, t) ■ d aP+ (a, t)) = j ^ 2 H{k, t) = d x (G^r 2 ^ 1+2 ^t- 2 , (3.7) 

5 In the absence of gravitational wave smoothing, these could even be at the scale set by the thickness of 
the string [3]. 
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where 

d x = Ac x l -^T[2{1 - X )(l + 2 X )\ [b x (l + x)(l + 2 X )] 2 ^ 1 )( 1+2 ^ (3.8) 

takes the value 0.006 in the matter era and 0.008 in the radiation era. 

On the long strings in the network, if the functions p+(u) and p~(v) were smooth and 
varied on the scale of the correlation length (which is not far below the horizon scale), then 
the typical velocity of any crossing would be the inverse correlation length, and this would be 



the size of cusps. What we see from Eq. (3.8) is that such slow crossings and large cusps are 
impossible. Superimposed on the slow motion is the irregular motion from the short distance 
structure, so all crossings occur at the much larger velocity V ~ 0.1(G/i)~ ( ^ 1 ~ x - )( ' 1+2x - ) t~ 1 , and 
all cusps have a size of order 

V- 1 ~ 10(G» (1 - x)(1+2x) t , (3.9) 

which is very much smaller than the correlation length. Effectively each large cusp breaks 
up into a large number of small cusps: the high fractal dimension of p+(u) and p_(f ) above 
the gravitational radiation length, noted in Ref . [8] , implies that these curves will cross many 
times near the would-be large cusp. We can estimate the number: the total length of each 
curve p + (u),p-(v) during one Hubble time is Vt ~ 0.1(G7/)- (1 - x)(1+2x) . The total number 
of crossings is roughly the total length of either path divided by the solid angle 4ir of the 
sphere, and so up to numerical factors it is large to the same extent that the size of each 
cusp is small. We could also imagine course graining, looking at the string with a resolution 
5a that is longer than the gravitational radiation scale, and the same argument would imply 
0((t/5<T) 1_x ) cusps each of size 0{t x 8a x ~ x ). 

Recall further that the cusp region is the site of rapid loop production, so these small 
cusps will actually end up on small loops, disconnected from the long string. It might seem 



paradoxical that the cusp size (3.9) is parametrically larger than the loop size (3.1), but this 
is another Lorentz effect: in the loop rest frame these sizes are of the same order. 

Regarding the large loops, we have already noted that these are likely to form without 
cusps. However, when they eventually decay by gravitational radiation at t ~ 1/TGfM, the 
higher harmonics will decay more rapidly, smoothing out the string and likely leading to new 
cusps. Thus there is the possibility of a cusp signal [221 [23] from these, although it will be 
reduced to the extent that the large loops represent only a fraction of the total loop length. 
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